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(57) ABSTRACT 

The invention is a method for obtaining the displacement of a 
flexible structure by using strain measurements obtained by 
stain sensors. By obtaining the displacement of structures in 
this manner, one may construct the deformed shape of the 
structure and display said deformed shape in real-time, 
enabling active control of the structure shape if desired. 
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Non-uniform cantilever beam instrumented with bending and distortion strain sensors. 
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METHOD FOR REAL-TIME STRUCTURE 
SHAPE-SENSING 

STATEMENT OF GOVERNMENT INTEREST 

5 

The invention described herein may be manufactured and 
used by or for the Government of the United States of 
America for governmental purposes without payment of any 
royalties thereon or therefor. 

10 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

This invention relates in general to determining the 
deformed shape of complex structures using stain sensor 15 
measurements in real-time and more particularly to employ- 
ing closed-form analytical solutions in conjunction with 
strain sensors to determine the deflection of complex struc- 
tures of real-time. 

2 . Description of the Related Art 20 

Resistive strain gauges have been employed for many years 

in experimental stress analysis in order to determine the 
stress-state produced in structures subjected to thermal and 
mechanical loads. More recently other sensors, such as those 
employing fiber optics, including fiber Bragg grating sensors, 25 
have also been used for the same purpose and for structural 
health and usage monitoring. Strain sensors measure the local 
strain in the plane of the structural surface to which the 
sensors are attached. These types of sensors, in conjunction 
with simple algorithms provide a straight-forward and trac- 30 
table methodology for determining the stress in structures. 
Conversely, a simple and straight forward approach for deter- 
mining the out-of-plane structural displacement is not avail- 
able, especially for lightweight highly flexible and complex 
structures. 35 

Therefore, it is desired to provide a method of obtaining 
both in-plane strain and out-of-plane displacements simulta- 
neously for complex, flexible structures. Further, it is desired 
to obtain these measurements in real-time. 

40 

SUMMARY OF THE INVENTION 

The invention proposed herein comprises a method for 
obtaining the out-of-plane deformations (i.e. structural defor- 
mations due to bending) of a flexible structure by using strain 45 
measurements obtained by strain sensors. By obtaining the 
displacement of structures in this manner, one may construct 
the deformed shape of the structure and display said 
deformed shape in real-time, enabling active control of the 
structure shape if desired. 50 

Accordingly, it is an object of this invention to provide a 
method for obtaining displacement of a flexible structure. 

It is a further object of this invention to provide a method of 
using strain sensor measurements to obtain displacement of a 
flexible structure. 55 

It is yet a further object of this invention to provide a 
method to obtain near real-time displacement of a flexible 
structure. 

This invention meets these and other objectives related to 
deformation measurements by providing a method for obtain- 60 
ing displacement of a structure having an axis and length. The 
method comprises the steps of first, dividing the structure into 
a plurality of sections having equal lengths. Strain sensors are 
placed adjacent to each section so one can get strain measure- 
ments related to each section of the structure. The strain 65 
measurements for each section of the structure are used to 
calculate the displacement by using analytical closed-form 


2 

methods (classical bending beam theory) that employ the 
length of the structure, the distance of the strain sensors from 
the axis of the structure, and the section length. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the drawings, which are not necessarily to scale, like or 
corresponding parts are denoted by like or corresponding 
reference numerals. 

FIG. 1 depicts a cylinder structure to which the present 
invention method may be employed. 

FIG. 2 depicts a wing-box structure to which the present 
invention may be employed. 

FIG. 3 depicts an aircraft fuselage structure to which the 
present invention may be employed. 

FIG. 4 depicts the results of testing the present invention on 
a cantilever beam structure. 

FIG. 5 depicts a non-uniform cantilever beam structure to 
which the present invention may be employed. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

The invention, as embodied herein, comprises a method of 
obtaining displacement of a flexible structure by using strain 
measurements taken along the structure. For the present 
invention, the term displacement is defined as deformation of 
a structure out of the plane of the structure as applied to the 
structure’s neutral axis. The term neutral axis is defined as the 
axis obtained by determining the axis or center point at each 
cross section of the structure. 

Many technologies that employ flexible structures, particu- 
larly those that are relatively light-weight, can benefit from 
near real-time displacement data. These include high-aspect 
ratio flexible structures, aircraft wings (e.g. UAVs), rotocraft 
vehicles, space vehicles, wind turbine blades for the alterna- 
tive eneigy community, space-based structures (booms and 
antenna), and long-span civil structures (bridges, dams). 
Since local strain is used as an input to the structural defor- 
mation method described herein, this method serves multiple 
purposes of structural shape sensing, loads monitoring, and 
structural health monitoring. 

Displacement Theory 

In general, the present invention employs a method of 
converting strain data obtained via strain sensors to get dis- 
placement data for flexible structures. To set forth the general 
theory of the conversion, the basic displacement equations 
will be developed for the weak nonuniform cantilever beam 
(e.g. Helios wing uniform tubular spar as limit case) subjected 
to bending, torsion, and combined bending and torsion are 
developed and described below. 

Bending 

The formulation of the displacement theories is built upon 
the classical beam differential equation (elastic curvature of 
deformed uniform beam) is given below 


tfy _ Mix) (D 

dx* ~ El 


where y is the vertical displacement, x is the span-wise coor- 
dinate, M(x) is the bending moment, E is the Young’s modu- 
lus, I is the moment of inertia. 
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At any beam cross-section x, t be bending moment M(x) in 
equation (1) may be related to the bending strain €(x) at the 
beam bottom (or top) fiber as 


tan*?,- = r 

1 


fW 

c(x) 


dx -f tan0,_i 


(7) 


€{X) 


M{x) s £/ 


where c is the beam half depth. 

In light of equation (2), equation (1) becomes 


fW 

rf* 2 C 


( 2 ) 


in which tan 0^ is the beam slope at the sensing station x f _ 1 . 

Substituting the beam depth equation (5) and strain equa- 
tion (6) into the slope equation (7), and after performing 
10 integration, one obtains 


(3) 


15 


tanei= 2^r[( 2_ ^) ei - i+e '] +taiiei -‘ 

(/ = 1, 2, 3, ... ,«) 


( 8 ) 


Note that, under the present strain formulation, the beam 
differential equation (3) contains only the beam half depth c 
and the bending strain e(x). The flexible rigidity term El is 
eliminated. The beam differential equation (3) for the uni- 
form beam could be used with sufficient accuracy for weak 
non-uniform beams. Namely, the beam half-depth c is no 
longer constant but a weak function of x [i.e., c==c(x)]. Hence, 
for the weak non-uniform beam, equation (3) may be modi- 
fied to 


20 


When i=l, tan 0 o (=O) will be the slope at the built-in end 
reference sensing station Xq. 

Deflection Equations 

The deflection y, at the strain sensing station x, may be 
obtained by carrying out double integration of equation (4) as 


25 




(9) 


£vy = fW 

dx 2 c(x) 


(4) 


30 


35 


The formulation of the present displacement theory for the 
non-uniform beams is based upon the modified beam differ- 
ential equation (4). 

The non-uniform beam may be considered as a piecewise 
tapered beam (either tapering down or tapering up within 
each segment). FIG. 1 shows the uniform cantilever beam 
(special case of non-uniform beam) with length 1 and weakly 40 
varying beam-depth 2c(x). The beam is installed with n+1 
equally spaced strain sensors at the sensing cross station x* 
(i=0, 1 , 2, 3, . . . , n) on the bottom surface of the beam. Thus, 
the beam is divided into n sections of equal length, Al”l/n. 
Each beam section may then be considered as a linearly 
tapered beam section (either tapering down or tapering up 
within the section). Therefore, within the region x / _ i <x<x < 
between the two adjacent sensing stations {x,_ l5 x,} at Al(=x,- 
x^l/n) distance apart, the beam half-depth c(x) and the 
bending strain e(x) may be expressed as linearly decreasing 
(or increasing) function of x as 


45 


in which y fml is the beam deflection at the sensing station x / . 1 . 

In light of the beam depth equation (5) and strain equation 
(6), the deflection equation (9) may be integrated twice to 
yield 


(Alf JY C ; \ -| 

» = 6^ll 3 ■ H 4 -*- 1 + 

(i = 1,2,3 ...«) 


( 10 ) 


When i=l, y 0 =tan 0 o =O at the built-in end reference sens- 
ing station x 0 . 

Combining equations (8) and (10), and using the recursion 
relationships, the deflection equation (10) becomes: 

Non-Uniform Beams 


(ID 


50 


^Zi{[ 3(2y - i) - (3y_2) ^K 4+w " 2 ^ , } H 


c{x) = C/-1 - (C/_1 - Ci) Xi~ t < X < Xi 


c{x) - - (©/_! - Si ) — — — ; x { -i <x<Xi 




(5) 


(6) 


>o 4- / A/tan^o 


55 


0 = U2,3 n) 


Uniform Beams (Limit Case) 


60 

where {c M> c,} and e,} are respectively the values of 
c(x) and e(x) at the sensing stations {x,_ l3 x,}. 

Slope Equations 

The slope tan 0, at the strain sensing station x, may be 65 
obtained by integrating equation (4) once between x^ and x, 


y, = 


(A0 2 


6c 


(/=1,2, 3 n) 


+ yo + tan^o 


( 12 ) 
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Torsion 

Similarly for torsion, Let T* be the twisting moment at the 
distortion sensing station x,, and let y, be the associated local yi 

surface angle of twist (shear strain) at the distortion sensing 
station x, in relation to the distortion sensing station x *„ t (FIG. 5 ^2 

3), then the torque/distortion relationship may be written as ^ “ "6^“H 3 ~ 




c n 


6 C0 lV 3_ ^ + £l J^° +A/tan ^ 


(17- a) 
(17-b) 


T)c 

GJ 


(13) 


10 


m \u 

yi= &rl( 3 - 


where G is the shear modulus, and J is the polar moment of 
inertia, and c is the outer radius of the structural member. The 15 
surface distortion y t (shear strain) is to be obtained from the y4 _ j£j > ffa 
distortion sensors as 603 

(Ai) 2 


[(l^ “ 7“) fi0 + ] + > 0 + 3A/tan0 o 


(17- d) 


6cz 


9- 


« = %^=zH= 2 < 1+ ^ 


(14) 


20 


_ (Afl 2 r/ 

The cross sectional twist angle ^ at sensing cross section x* " 6c 4 l\ “ 

may then be calculated from 25 (a/) 2 ,, 

6c 2 n 


M 

<f>; = — (70 + n + 7z + ■ ■ ■ + r/-i) 


(15) 


(17 ' c) 
6cq 

4 Sh*^] + S ls - 7 ?) , »* 7 *] + 

^ [(21 - lO^jeo + 10c,] + jo + 4&haaBa 

(17 ‘ e) 

[(27 - 13^-)^> + 13fi! ] + > 0 + 5Man$, 


30 


Combined Bending and Torsion 
In order to obtain the true bending strain the value of 
measured strain €* must be corrected by using the following 
bending strain correction equation. 35 


At the beam - tip (i = n ): 

m 2 


y« = 


6c 


7=1 


3(27-1)- (3/ -2) 


c n-J+l 




(17- f) 


(3 j - 2)e n _ J * + i \ + yo + ttAftau0 o 




COS0/COS7/ 


(16) 


where the surface twist angle y t and the cross sectional twist 
angle (j),. are to be calculated respectively from previous equa- 
tions and input into the following equation 

Non-Uniform Beams 


;-2)z lW+ ,j 


(17) 


(3y- 

1,2,3,... ,*) 


+ yo + i’A/tan^o 


Uniform Beams (Limit Case) 


. (A if 

6c 


(3/-1)£o+-6^ {i~jySj+?i 
7=1 


+ jo + tan9o 


(18) 


Writing out equation (17) explicitly for different indices i 
for displaying functional behavior, there results. 


40 in which the built-in condition y 0 =tan 0 o =O holds. 

Equations (17-a)-(17-f) explicitly show that the deflection 
y t (i=l , 2, 3, . . . , n) is obtained by summing up the geometrical 
and strain data evaluated at the inboard strain sensing loca- 
tions (x 0 , x x , x 2 , . . . , x,) including the current strain sensing 
5 location x,. By obtaining the deflection at all of the strain 
sensing locations, one may thereby generate the deformed 
shape of the tapered cantilever beam such as an aircraft wing. 

It is important to mention that the deflection equation (17) 
5Q developed for slightly non-uniform beam can now be used for 
the uniform beam case by setting c^j^c, as a limit case. 

In general, the present invention includes the following 
steps: 

1 . A structure is selected in which strain and displacement 
55 measurements are required. These can be classic structural 
configurations, such as beams, plates, shells, tubes, etc. or 
structures of complex configurations, such as tapered beams, 
wing foils, space-based antenna, wind turbine blades, 
rotocraft blades, and other complex structures. 

60 2 . Discretize the spatial domain of the structure by dividing 

the structure into sections, spaced at constant distances apart. 

3. Install strain sensors on the surface of the structure so 
that the sensors are located adjacent to each of the sections 
defined in step 2. The strain sensors preferably are placed 

65 equidistant from each other. 

4. Measure the strain at each section using the strain sen- 
sors located along the surface of the structure. 
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5. Convert the strain data into displacement data employing 

analytical closed-form methods (classical beam bending 
theory) that employ the length of the structure, the distance of 
each strain sensor from the neutral axis of the structure at its 
location, and the section length. 5 

6. Optionally plot the deflection data at the strain sensor 

locations to construct the deformed shape of the structure to 
display the deformed shape in real-time or to input the 
deformed shape into a control system to actively control the 
shape of the structure. 10 

In order to practice the present invention, strain sensors 
must be placed on the structure at issue. While any type strain 
sensors that accurately measure strain on a structure (for 
example, strain gauges), may be employed in the present 
invention, it is preferred that a light-weight strain sensor, such 15 
as those using fiber-optic technologies, be used. This is 
because many of the structures named above, that might 
employ the present invention, are relatively light-weight and 
flexible themselves. 

Preferably, fiber optic sensors are employed in the present 20 
invention. The most preferred strain sensor to practice the 
present invention are Fiber Bragg Grating (FBG) sensors. 
These sensors are preferred because they are minimally 
obtrusive, ultra-lightweight, easily installed, accurate, 
immune to EMI, and inherently-safe (no joule heating, spark- 25 
ing). The most preferred configuration and use of FBG sen- 
sors for the present invention employ the Optical Frequency 
Domain Reflectometry (OFDR) technique with hardware 
architecture described in U.S. Pat. No. 5,798,521 and U.S. 

Pat. No. 6,566,648 which are incorporated herein by refer- 30 
ence. This approach uses low reflectivity gratings all with the 
same center wavelength and a tunable laser source. The FBGs 
are preferably located on a single optical fiber. This allows 
hundreds of strain sensors to be located down the length of the 
fiber, a common configuration is to use 480 FBGs on a single 35 
fiber spaced at 1 FBG/cm. This configuration allows strain 
measurements to be acquired at much higher spatial resolu- 
tion than other current senor technologies, making it flexible 
enough to employ a user-selected grating density depending 
on the type of application. 40 

The general approach described herein may be employed 
on both simple and realistic, complex structures, such as 
aircraft wings, as well as a wide-variety of complex structures 
which vary spatially in all three Cartesian coordinates (x,y,z), 
for example tapered beams, which vary in all three axes, 45 
plates of nonuniform thickness, hollow tubes, beams, col- 
umns, shells, etc. 

Also, this approach, unlike many prior art approached for 
gauging structural changes, is not material independent. It 
does not depend of the material properties of the substrate 50 
such as modulus of elasticity. Poisson’s ratio, bulk modulus, 
etc. Therefore no pre-test calibrations or advanced knowledge 
about these complex conditions is needed to practice the 
present invention. 

The method of the present invention is also completely 55 
independent of applied loads (concentrated, distributed) and 
boundary conditions being applied to the structures at issue. 
The invention can operate on structures under any general 
loading conditions including moments/forces, which are spa- 
tially nonuniform. The invention may also be used on any of 60 
classical beam/plate/shell boundary conditions such as fixed 
ends, simply-supported, flexible supports, etc, or any general 
condition in-between without any dependence on heteroge- 
neous boundary conditions in order to simplify the closed 
form solution. 65 

Further, the invention does not apply the typical assump- 
tions required to solve classical analytical linear partial dif- 
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ferential equations. Such assumptions significantly reduce 
the types of structures that may be assessed. 

A specific example of the present invention is described 
below. Referring to FIG. 1, using a simple cylinder 100 of 
length, 1, and radius, c, the sections 102 are defined to be A1 
distances apart. 

Strain sensors 104 (e.g. strain gages, fiber Bragg gratings, 
etc.) on the surface of distance c from the neutral axis 106 
down the length of the structure 100 so that a strain sensors 
104 are located adjacent to the sections 102. 

Assuming a specific case in which the cylinder 100 is 
subjected to lateral loads only, which produce bending in the 
cylinder 100, then the displacements due to bending at y f can 
be calculated by inputting the cylinder’s 100 length, 1, the 
distance the sensors 104 are away from the neutral axis 106, 
c, the spacing distance A3, and the measured strains at each 
section 102, €„ into the following generalized equation in 
which the thickness of the structure varies from station to 
station in the x direction (i.e. c is nonuniform as a function of 

x) 

y- t = (*9) 

¥Zi{[ 3w - i) - (3/ - 2) ^ i l e ^ +w - 2 ^‘} + 

y=l 

yo + /A/tan0o 

1,2, 3, n) 

Since c(x)=constant for cylinder 100, then equation 19 
reduces to 


(A/) 2 £4 

yi ~ (3i-l)£o + 6^ (i-fej + Ci +■ yo + tanflo 

M 

(Is 1,2, 3, ... ,«) 


By plotting the deflections y„ (i-1 , 2, 3, . . . , n) at the strain 
sensing sections x„ the deformed shape of the cylinder 100 
may be constructed, and could be either be displayed for real 
time deformation status or input into the control system to 
actively control shape. 

The above example assumes that the out-of-plane displace- 
ments being determined are produced by loads applied nor- 
mal to the structural surface and that the in-plane strains are 
very low. This assumption is reasonable for a plethora of 
structural configurations such as wings, cantilevered beams, 
turbine blades, pressure bulkheads, etc. If this assumption is 
not valid, the invention may still be used by employing two 
sensor sets located on the top and bottom surfaces of the 
structure. This configuration allows a user to take the mean of 
the strain measurement of the top and bottom sensors at each 
section of the structure and input the mean value into the 
appropriate conversion formula. This technique then effec- 
tively eliminates the unwanted in-plane strain contribution at 
each section number produced in the structure by in-plane 
loads. 

FIG. 5 shows a non-uniform cantilevered beam. Equation 
1 9 would be employed to practice the present invention in lieu 
of equation 20, as described in the above example. 

In addition to the examples of conversion equations devel- 
oped and discussed above, the deformations of a multitude of 
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other structures can be solved using the same basic equations 
19 or 20 above. Following are examples for complex struc- 
tures. 

APPLICATION TO WING BOXES 

For the installations of strain sensors for a wing box-like 
structure under combined bending and torsion, two bending 
strain sensor lines may be used each as input to the basic 
equation. FIG. 2 shows a wing box under combined bending 
and torsion. The wind box is a type of a slightly tapered 
cantilever beam. The two bending strain sensing lines system 
is a simple way to simultaneously monitor the wing box 
deflections and the wing box cross sectional twists. The two 
(front and rear) bending strain sensing lines 310, 312 are 
parallel, separated by a chord-wise distance of d, and are 
oriented in the span-wise direction. 

1 . Deflections 

The deflections {y*, y’,} of the front and rear sensing sta- 
tions 310, 312 lying at the same span-wise sensing cross 
section \ (i=0, 1, 2, 3, . . . , n) may be calculated from 
deflection equation (19) which is rewritten in the following 
for the front and the rear strain sensing stations. 

For the front strain sensing line 310 equation x is: 

yi = (21) 

b i 

>o + (OAltan^o 

(i = 1, 2, 3, ... , n); yo = Man0o - 0 
For the rear strain sensing line 312 equation x is: 
yi = (22) 

;=t 

y'o + (i)A/tan0o 

(i= 1, 2, 3, ... ,n); / o = Man^ o = 0 


2. Cross Sectional Twists 

If d denotes the chord-wise distance between the front and 
rear bending strain sensing lines 310, 312, then the cross 
sectional twist angle ^ at any strain sensing cross section x, 
may be calculated from the following cross sectional twist 
angle equation. 

<Pi = (i = 1, 2, 3, ... ,») (23) 


in which <|> 0 =0 at the wing root. The deflections {y,, y'J may 
be calculated respectively from the deflection equations (21), 
(22) using the measured bending strain data as inputs. With 
the use of equation (23), installations of the distortion sensors 


10 

are not required. An alterative method for the predications of 
the cross sectional twist angle <j>, for the wing box is discussed 
in the following example. 

5 APPLICATION TO AIRCRAFT FUSELAGES 

Referring to FIG. 3, an aircraft fuselage is a tubular shaped 
structure with varying cross sections. This kind of structure 
may be considered as a non-uniform unsupported free-free 
beam during flight. To measure the deformed shape of the 
airborne fuselage due to vertical and horizontal bending, and 
torsion, three sensing lines 420, 422, 424 are required. The 
bottom 420(belly) and side 422 sensing lines are to monitor 
15 the vertical and horizontal deflections respectively. The tor- 
sion sensing line (not depicted) may be installed either on the 
bottom (or side) surface parallel to the bottom (or side) bend- 
ing sensing line. 

In formulating the deflection equation (19) for a slightly 
20 non-uniform cantilever beam, duplicated below as equation 
(24): 


x = < 24 > 

^E^r{[ w - 1) - <3/ - 2) ^ L b +(3; - 2)e ^ , } + 

h i 

y 0 + (i)Man0 o 

(/ s= 1 , 2, 3 n) 

the reference strain sensing station x 0 is located at the 
35 built-in end. With the built-in condition y 0 =tan 0 o =O imposed, 
equation (24) is also applicable to the moving cantilever 
beams because, the deflection y t calculated from equation 
(24) is relative to the tangent line (tan 0 o =O) stemming from 
the reference strain sensing station Xq at the built-in end. 
40 Equation (24) may be applicable to the free-free unsupported 
moving fuselage case in the following fashion. 

The fuselage may be considered as a free-free unsupported 
beam consisted of two non-uniform cantilever beams joined 
together at the center-of-gravity (CG) cross section. During 
45 flight, the CG cross section has the least movement compar- 
ing with other fuselage cross-sections. For practical purpose, 
it is ideal to choose the reference strain sensing station Xq at 
the CG cross-section, and divide the strain sensing line into 
two segments 430, 432. One segment is for the forward- 
50 fuselage section and the other for the aft-fuselage section. 
Thus, by imposing y 0 =tan 0 o =O at the reference sensing sta- 
tion Xq (considered as a built-in end), deflection equation (24) 
for the cantilever beam may be used to calculate the deflec- 
tions of the fore- and aft-fuselage sections. The calculate 
55 deflection y t calculated from the deflection will then be the 
relative deflection with respect to the tangent line passing 
though the reference sensing station x 0 , which moves (trans- 
lates and rotates) with the fuselage. 

The present invention was tested on a 10-fit long, thin- 
60 walled composite beam. Results are shown in FIG. 4 from 
cantilever beam testing for three combined load cases: 2-lbs 
bending/10 in-lbs torque, 4-lbs bending/20 in-lbs torque, and 
6-lbs bending/20 in-lbs torque. For each load case, the deflec- 
tion of the cantilever beam is shown using 1) strain gage 
65 measurements input into the strain/displacement equations, 
FBG strain sensors input into the strain/displacement equa- 
tions, and direct measurements using potentiometers. Good 
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agreement is shown in FIG. 4 for all three cases, thus validat- 
ing the equations for cantilever beam deflection developed 
above. 

What is described are specific examples of many possible 
variations on the same invention and are not intended in a 5 
limiting sense. The claimed invention can be practiced using 
other variations not specifically described above. 

What is claimed is: 

1. A method for obtaining displacement of a structure 
having a depth with a neutral axis and length, comprising the 10 
steps of: 

dividing the structure into a plurality of sections, the sec- 
tions having equal lengths; 

providing strain sensors adjacent to each section; 

obtaining strain measurements horn each strain sensor; 15 
and, 

calculating the displacements on the structure out of a 
plane from the neutral axis using the strain measure- 
ments, the structure length, the strain senor distance 
from the structure neutral axis, and the section length. 20 

2. The method of claim 1, wherein the strain sensors com- 
prise a location at a first end of each section, the locations 
defined as Xq, x 15 x 2 . . . x„, wherein n is the number of strain 
sensors and x 0 is located at a first end of the structure, x„ is 
located at a second end of the structure, and the remaining 25 
strain sensors are located in numerical order between x 0 and 
x„. 

3. The method of claim 1, wherein the strain sensors com- 
prise fiber optic sensors using a tunable light source. 

4. The method of claim 3, wherein the strain sensors com- 30 
prise fiber Bragg grating sensors. 

5. The method of claim 4, wherein the fiber Bragg grating 
sensors are located along a single optical fiber on a surface of 
the structure. 

6. The method of claim 5, wherein the sections comprise a 35 
minimum length of about once centimeter. 

7. The method of claim 1, wherein the providing strain 
sensors step includes providing strain sensors adjacent to 
each section at a top and bottom location. 

8. The method of claim 7, wherein the strain sensors are 40 
located along two optical fibers, a first optical fiber on a top 
surface of the structure and a second optical fiber on a bottom 
surface of the structure. 

9. A method for obtaining displacement of a structure 
having a depth with a neutral axis and length, comprising the 45 
steps of: 

dividing the structure into a plurality of sections, the sec- 
tions having equal lengths; 

providing strain sensors comprising a location at a first end 
of each section, the locations defined as x 0 , x l5 x 2 . . . x„, so 
wherein n is the number of strain sensors and x 0 is 
located at a first end of the structure, x„ is located at a 
second end of the structure, and the remaining strain 
sensors are located in numerical order between Xq and 


obtaining strain measurements from each strain sensor; 
and, 

calculating the displacements on the structure using the 
strain measurements, the structure length, the strain sen- 
sor distance from the structure neutral axis, and the 
section length further comprising calculating deflection 
at each strain sensor, defined as location x*, by summing 
the strain measurement at location x* with the strain 
measurement at strain sensors at locations x 0 through 
x,-i* 

10. The method of claim 9, wherein the calculating step 
comprises the use of the algorithm: 


yo + 1 A/tando 


where y t is the displacement at a strain sensor locationx*, A1 
is the distance between two adjacent strain sensors, 1 is 
the length of the structure, c is the distance from the 
strain sensor location and the neutral axis, and € is the 
measured strain at each strain sensor location. 

11. The method of claim 10, wherein the distance from 
each strain sensor to the neutral axis is equal, resulting in the 
algorithm: 


W 
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12. The method off claim 10, wherein the strain sensors 
comprise fiber optic sensors using a tunable light source. 

13. The method of claim 12, wherein the strain sensors 
comprise fiber Bragg grating sensors. 

14. The method of claim 13, wherein the fiber Bragg grat- 
ing sensors are located along a single optical fiber on a surface 
of the structure. 

15. The method of claim 14, wherein the sections comprise 
a minimum length of about once centimeter. 

16. The method of claim 10, wherein the providing strain 
sensors step includes providing strain sensors adjacent to 
each section at a top and bottom location. 

17. The method of claim 16, wherein the strain sensors are 
located along two optical fibers, a first optical fiber on a top 
surface of the structure and a second optical fiber on a bottom 
surface of the structure. 
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